Transport through quantum dots with magnetic impurities 
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We analyze the electronic transport through a quantum dot that contains a magnetic impurity. 
The coherent transport of electrons is governed by the quantum confinement inside the dot, but is 
also influenced by the exchange interaction with the impurity. The interplay between the two gives 
raise to the singlet-triplet splitting of the energy levels available for the tunneling electron. In this 
paper, we focus on the charge fluctuations and, more precisely, the height of the conductance peaks. 
We show that the conductance peaks corresponding to the triplet levels are three times higher than 
those corresponding to singlet levels, if electronic correlations are neglected (for non-interacting 
dots, when an exact solution can be obtained). Next, we consider the Coulomb repulsion and the 
many-body correlations. In this case, the singlet /triplet peak height ratio has a complex behavior. 
Usually the highest peak corresponds to the state that is lowest in energy (ground state) , regardless 
if it is singlet or triplet. In the end, we get an insight on the Kondo regime for such a system, and 
show the formation of three Kondo peaks. We use the equation of motion method with appropriate 
decoupling. 
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I. INTRODUCTION 

The quantum dots are zero-dimensional structures (so called artificial atoms) that arc possible to be experimentally 
obtained and investigated individually, for some years now. An incoming electron may tunnel through a quantum dot 
only if its energy matches a dot energy level, so we have a single electron transistor [ij . Another interesting potential 
application is quantum computing using the spin of the quantum dot as information carrier . The parameters that 
influence transmittance through a quantum dot (like the size of the dot or the hybridization with the leads) are much 
more easy to control comparing to the parameters of a bulk system, so one can test quantum theories in a way that 
was not possible before. We give as example the observation of the two-channel Kondo effect, described theoretically 
many years ago 0, 0] , but realized experimentally only very recently, in a double-dot system Q . 

In the last years, people have started to address the problem of quantum dots with magnetic impurities 0, 0, H, S 
[13, [ll|, , which is also the subject of the present paper. The transient electron feels the exchange interaction with 
the impurity and forms singlet /triplet entangled states before tunneling out of the dot. We analyze the way in which 
this physics is seen in transport. It will be shown that the peaks corresponding to the triplet levels are three times 
higher than those corresponding to singlet levels, if we neglect the Coulomb repulsion and electronic correlations in 
the dot. In large dots, with low confinement, the Coulomb repulsion is reduced, and the approximation to neglect it 
may still capture some important physics. Such a case will be addressed in Section II. 

The case of small quantum dots, with strong Coulomb repulsion is treated in Section III. If the Coulomb interaction 
and the many-body correlations are considered, the picture becomes more complex: the peak corresponding to the 
ground state is usually higher, regardless if it is singlet or triplet. The peaks height depends on the ratio J/F (exchange 
interaction versus coupling strength with the leads). An analytical formula is proposed for the transmittance and in 
particular for the singlet/triplet peak height ratio. In the end, we get an insight on the Kondo regime and show the 
- formation of three Kondo peaks in the density of states. 
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II. LARGE QUANTUM DOTS WITH MAGNETIC IMPURITIES 

In this section we shall neglect all electron-electron interactions, except for the interaction with the magnetic 
impurity. The correlations are neglected as well, and they will be considered in the next section. The problem 
presented in this section is a two-electron scattering problem (one of which -the impurity- is fixed and may only 
change its spin orientation). 

If correlations are neglected, one can easily afford to consider a many-site dot, which is realistic from geometrical 
point of view. In (large) quantum dots with low confinement, the Coulomb repulsion is reduced, and this justifies the 
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non-interacting model used here. A general lattice Hamiltonian, with a localized spin interaction, can be written as: 



+ 



\j (cl^c„| - c\^Cni)S^^ + i J (4 ^C„|5„ + c\^Cn^S+) ] n & QD , (1) 



where c| (q) are creation (annihilation) operators in the dot sites indexed by i ; the index n is devoted to the site of 
the dot where the magnetic impurity is placed. Wij,o- are hopping parameters (actually the hopping parameters are 
spin-independent and in the following we give up the spin index) and are the phases associated with an (eventually 
applied) magnetic field. The last two term account for the interaction with the magnetic impurity [T^. 

From this point on, we take S— 1/2. One can define the fermion operators {dj, d|, d^, dj} for the localized spin 
(provided one projects out all the states with occupancy different from 1) 

Then, S^, S^and S~ in eq (1) can be written as: 

S:, = ^{d\d^~ ~ d\d^), S+=d\d^, S-=dld^. (2) 

The natural formulation of the problem is in terms of singlet- triplet operators (see for instance [13 )■ Let us introduce 
the singlet operator and the triplet operators T[ (p=l,2,3): 

E; = {d^Cii - dici^), Tl = (d|C,x +d4C,t), 

Tf = d|C,-f, Tf = d^c^i . (3) 

The above operators may be used in order to write the Hamiltonian (Eq.l) in the following way (where we have used 
also Ud = + = 1): 
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= e''^'^*'^- + S^jVg - - J(5.j5„, 



'''^*'^+^..K, + 7J%<S.„ . (4) 
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In order to analyze transport properties, one needs to connect leads to the quantum dot. The leads are modeled by a 
one-dimensional chain , consistent with the dot tight-binding model. For the non-interacting dot (more precisely: the 
interaction is only with the magnetic impurity, and the electron-electron interactions and correlations are neglected), 
the introduction of leads is equivalent to the introduction of a complex selfenergy in the Hamiltonian, at the contact 
sites. The procedure is described in detail in 0. It results the effective Hamiltonian: 

= E [(^.g+5.a-5,or^e-'^)stE, + (z«5+5,„5,„r^e-'^) ^ T^T^] 

= Hl/f + YH'T'/- (5) 

We describe briefly the next steps. First, we employ the Landauer-Buttiker formula to connect the transmission 
through the dot with the retarded Green function: 

Ta,a'S,S'{E) = 4T^stn\ I < a,S\G+^,{E)\a',S' > \\ (6) 

The above formula assumes equal coupling strength to the leads (tq) connected with the sites (states) a, a'. "S"' 
stands for the spin state of the magnetic impurity that is changed to " S"" after the the scattering . The electron spin 
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is also changed from a to a' and only the transitions that conserve the total spin give non-zero transmittances. The 
parameter "g" is the impulse of the incident electron (see and references therein for details about the dispersion 
formula in the leads and connection with the Landauer-Buttiker formalism) . 

At this point we remember that, for a single electron case, the retarded Green function is equal to the resolvent, 
and can be calculated by: 

G+{E)^{E-H'^-^f)-^. (7) 

The following definitions are introduced : 

Gt{E)^{E-H^^f)-\ 

GUE)^iE-H^V)-\ (8) 
All triplet Green functions are equal (we give up the index "p") and the transmittance can be written: 

Ttotai= J2 T^,.'S,S'iE)='ir^s^n\i\\GUE)\' + l\G+{E)\'). (9) 

a,cr'S,S' 

We show in Fig.l the transmittance peaks for a quantum dot modeled by a 3 x 5 two-dimensional lattice. The 
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FIG. 1: The conductance through a many-site quantum dot with a magnetic impurity, in the absence of electron-electron 
interaction (only the interaction with the magnetic impurity is considered) . Some peaks split into triplet and singlet, with 
the heights 3/4 and 1/4, respectively. The peaks corresponding to trajectories that avoid the impurity do not split and remain 
degenerate. 

exchange interaction causes the singlet-triplet splitting of the spectrum and the conductance peaks. The triplet peaks 
have the hight 3/4 and the singlet peaks 1/4. For a realistic 2D dot, certain cigenstates may be zero at the impurity 
position. Obviously, the corresponding transmittance peaks will not split. This case is also captured in Fig.l (the 
"degenerate peaks" do not undergo the singlet-triplet splitting). 

III. SMALL QUANTUM DOTS WITH MAGNETIC IMPURITIES 

We saw in the previous section that the non-interacting dot is easy to be modeled by a n x m lattice, but the problem 
becomes (technically) much more difficult when we include the Coulomb interaction. A general Hamiltonian, that 

includes the Coulomb interaction between any two electrons, can be written: : H = Hs + J2i j a a' crCi,crC^ a'^j,'^' ■ 

There are several ways to address such a complex Hamiltonian, that allows only approximate solutions with consider- 
able difficulty. In this section , we shall restrict to the case of a single-site dot with Coulomb interaction and exchange 
interaction with a magnetic impurity. In comparison to the Hamiltonian Eq.l, we shall introduce the Coulomb inter- 
action, but also we shall write explicitly, from the beginning, the Hamiltonian of the leads and the leads-dot hopping 
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term. The reason is the following: in the previous section we used the effective Hamiltonian trick, meaning that 
leads are introduced by a selfenergy term in the coupling sites. But this is no longer possible if we consider electronic 
correlations or the Kondo effect, when the role of the leads is more complex. 
We re-write the Hamiltonian, in a convenient notation: 

k,a.a. (7 

\ J (cJt^OT - cJiCOJ,)S"' + i J {cl^CoiS~ + cI^CotS+) + ^ ta{cl^Cka.a + h.c.) , (10) 

k,a.a 

where the first term represents the electrons in the lead a — L, R. The second term stands for the electronic level 
eo in the dot. The next four terms describe interactions: the Coulomb interaction of electrons with the opposite spin 
orientation at the level eo and exchange interactions with the magnetic impurity. The last term in the Hamiltonian 
(Eq.lO) corresponds to the coupling between the quantum dot and the leads. It shall be considered, for simplicity, 
th ~ tB. = t- 

Now, we want to determine the conductance for the model described by the Hamiltonian (Eq.lO). The current can 
be expressed by means of the non-equilibrium Green functions as [l^ [ig : 

J = ^^J dujihiuj) - /fl(^)](-Im((coTl4))), (11) 

where F = 2T:t'^p and p = 1/2D is the DOS for the square band approximation; the half-width D will be taken 
as unity. and fji are the Fermi distribution functions in the left and right leads, respectively. ((co||cJ|)) is the 
retarded single particle Green function for an electron with the spin a =t at the QD, which can be determined by 
the equation of motion (EOM). 




FIG. 2: The height of the singlet and triplet peaks, that depend on the ratio J/F (exchange interaction versus coupling strength 
with the leads) , if the Coulomb interaction and electronic correlations are considered. 

The generic equation of motion (EOM) for the energy dependent retarded Green function is given by 

u;{{A\B)) = {{A,B}) + {{[A,H]\B)) , (12) 

where {{A, B}) is the thermal average of the anticommutator between the operators A and B. The procedure for solving 
the EOM for the retarded green functions was described in detail in [T^j ■ Basically, the equations of motion introduce 
higher order Green functions, which result by performing the commutation (with the Hamiltonian) required by the 
formula above. For a finite system, even in the presence of interactions, one can write a finite number of independent 
equations and the system closes (even if sometimes the number of equations is very big, it is still solvable) . But for 
our case, interacting system -I- infinite leads, one has an infinite set of equations, and approximations are needed to 
close the system. 

We propose the next approximation (that is valid at high temperatures (iTj ) 

2tJ2{{c,^,A\clj) « -^^{{co,.A\clJ). (13) 

k 
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The following solution is obtained for the conductance through the singlet and triplet peaks (we take U oo) 



G = Gs + Qt 



1 + 202 



-R. 



3(1 + 205) 



2(3 + 40s - 40s0t) 2(3 + 40s - 40s0t) 



(14) 



where i?s = (2eV/i) x rVp^ + (Ep - eo + 3J/4)2] and Rt = (2eV/i) x T^/[r^ + {Ep - eo - J/4)2] corresponds to 
the resonant conductance through the singlet and the triplet level, respectively, 0s — arctan[(eo — Ep — 3J/4)/r]/7r 
and 0T = arctan[(eo — Ep + J/A)/T]/tt. The fractions that multiply the resonances i?s(T) contain the information 
about the height of the conductance peaks. The ground and excited states role changes with the sign of J, and the 
dependence on the coupling also changes, being much more pronounced for the excited state. All these effects are 
incorporated in the above formula. 

By employing further the simplified form Eq.(14), one can estimate the ratio between the heights of the conductance 
for the singlet and triplet peaks. We notice that in Eq.(14), the parameter 0s vanishes at the singlet resonance (found 
with the condition eo ~ Ep — 3 J/4 = 0). 0t vanishes at the triplet resonance. Therefore, one can straightforward 
calculate the ratio of the conductance peaks 



max[gs] _ (l + 20o)(3- 
ma.x[GT] 



9(1 - 20o) 



(15) 



where 0o = arctan(J/r)/7r. Fig. 2 gives a graphical image of the peaks height dependence on the parameter J/F. For 
the antiferromagnetic coupling the electronic transmission through the singlet state dominates and in the limit J 3> F 
this ratio goes to infinity and the triplet peak disappears. In the case of the ferromagnetic coupling and | J| ^ F we 
have opposite situation - the singlet states disappears from transport and only the peak corresponding to the triplet 
ground state is visible. 
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FIG. 3: The density of states for a quantum dot with a magnetic impurity, in the Kondo regime. Notice the formation of three 
Kondo peaks, one at the Fermi energy (Ep — 0) and two side peaks at a distance equal to the exchange interaction (J = 0.3) 



In the last part of the paper, we give a very brief analysis of the Kondo regime for the impurity dot. It is known 
that a Kondo peak in the density of states, at the Fermi energy, is present only if a degenerate level exists below the 
Fermy energy. This would allow for fluctuations (in the "classical" Kondo one has spin fluctuations), that involve also 
the electrons from the leads which, for this reason, can pass through the dot with increased probability. In our case, 
such a degenerate level would be the triplet and the system shows the central Kondo peak (at the Fermi energy). 
This central peak is important for transport, because at low temperature and low bias, practically only the electrons 
around the Fermy level contribute to the conductance. 

On the other hand, the singlet-triplet level structure should allow also the observation of side Kondo peaks, similar 
to the case of an applied magnetic field We obtain a three-peaks Kondo structure and the typical density of 

states is plotted in Fig. 3. The calculations were performed following the decoupling proposed in [l8l [l9| (see [l^l for 
calculations details) for temperatures closed to the Kondo temperature. 
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IV. CONCLUSIONS 



111 this paper, we present two models of quantum dots with magnetic impurities. The first model is a many-site 
lattice model (large quantum dots). It neglects the Coulomb interaction and many-body correlations, but can be 
used to model realistic geometries. In the absence of the magnetic impurity, the transmittance shows peaks with 
the height 1 (perfect transmittance). If a magnetic impurity is placed in the dot, one has triplet peaks with height 
3/4 and singlet peaks with height 1/4. Some eigenstates inside the dot may avoid the impurity position, and the 
corresponding peaks do not split into singlet and triplet. 

The second model we propose is a single-site dot with Coulomb and exchange interactions. Here we give up the 
geometrical aspects, in order to focus on the correlations effects. Such a model is realistic for small quantum dots, 
where the Coulomb is strong and the level-spacing is large. The singlet and triplet conductance peaks are shown to 
depend on the ratio J/T (the exchange interaction versus the coupling strength with the leads). In particular, the 
transmission through the excited states depends strongly on the coupling to the leads. 

In the Kondo regime, the presence of the 5 = 1/2 magnetic impurity generates three Kondo peaks in the density 
of states. 
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